Open charm production in relativistic nucleus-nucleus collisions * 
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Abstract 



I We calculate excitation functions for open charm mesons in Au+Au reactions 

I from AGS to RHIC energies within the HSD transport approach which is 

o ' 

O I based on string, quark, diquark {q, q, qq, qq) and hadronic degrees of freedom. 



The open charm cross sections from pN and vrA^ reactions are fitted to results 
from PYTHIA and scaled in magnitude to the available experimental data. 



O I From our dynamical calculations we find an approximate mj^-scaling for pions, 

^3 ' 



kaons, D-mesons and J/^ - when discarding final state elastic scattering of 
kaons and (/)-mesons with pions - in central collisions of Au+Au at 160 A-GeV 



X 

I (with an apparent slope of 176 MeV) without employing the assumption of a 

Quark-Gluon Plasma (QGP). We demonstrate that this result is essentially 
due to a relative mT^-scaling in pp collisions at ^/s ~ 17.3 GeV. At lower 
bombarding energies of 25 A-GeV a suppression of D-mesons by a factor of ~ 
10 relative to a global mr-scaling with slope 143 MeV is expected. However, 
when incorporating attractive D-meson self energies as suggested by QCD 
sum rules, an approximate mr-scaling is regained even at 25 A-GeV. The 
effects of D-meson rescattering and charmonium absorption are discussed, 
furthermore, with respect to rapidity and transverse mass distributions in 



*supported by GSI Darmstadt and FZ Jiilich 



central collisions of Au + Au at 25, 160 A-GeV and 21.3 A-TeV. 
PACS: 25.75.-q; 13.60.L2; 14.40.Lb; 14.65.Dw 
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I. INTRODUCTION 



Apart from the light and strange flavor {u, u, d, d, s, s) quark physics and their hadronic 
bound states in the vacuum (vr, K, etc.) the interest in hadronic states with charm flavors 
(c, c) has been rising continuously in line with the development of new experimental facilities 
This relates to the charm production cross section in pN and vrA^ reactions as well 
as to their interactions with baryons and mesons which determine their properties (spectral 
functions) in the hadronic medium. 

The charm quark degrees of freedom have gained vivid interest especially in the context 
of a phase transition to the quark-gluon plasma (QGP) [H where cc meson states should 
no longer be formed due to color screening However, the suppression of J/\l/ and 

mesons in the high density phase of nucleus-nucleus collisions P,^] might also be attributed to 
inelastic comover scattering (cf. |[T0|-[T5| and Refs. therein) provided that the corresponding 
J/\l'-hadron cross sections are in the order of a few mb [p!^-|2T|. Present theoretical estimates 
here differ by more than an order of magnitude [^] especially with respect to J/\E'-meson 
scattering such that the question of charmonium suppression is not yet settled. On the other 
hand, the enhancement of 'intermediate-mass dileptons' in Pb + Pb collisions at the SPS has 
been tentatively attributed to an enhancement of 'open charm' in nucleus-nucleus collisions 



relative to pA reactions at the same invariant energy ^/s . It should be mentioned that 
this enhancement does not stem from the charmonium dissociation since it is about two 
orders of magnitude larger the total charmonium yield. Thus 'charmonium suppression' and 
'open charm enhancement' are present facets of relativistic heavy-ion collisions. 

Furthermore, it is well known experimentally [|2^ that the D,D and D*,D* mesons 
show some analogy to the K,K and K*,K* mesons with respect to their excitation spec- 
trum because the strange (antistrange) quark is replaced by a charm (anticharm) quark in 
the hadronic state. Since quite substantial in-medium potentials have been suggested for 



antikaons in dense nuclear matter [p!0| , p5| , the latter might also show up for the correspond- 
ing D-mesons in view of a similar wavefunction for the light quark [^]. In fact, QCD sum 



rule studies point towards attractive potentials for the D-mesons which might lead to 
enhanced production cross sections of open charm especially at low bombarding energies 
close to threshold. Substantially lower in-medium effects are expected for the J/\l/ or t]c do 
to a small coupling of the c, c quarks to the nuclear medium . Thus the charm- meson sec- 



tor, which is insufficiently known so far, provides a theoretical |lT9| , |29| -p3[] and experimental 
challenge for the future ||3^j35[| . 

In this work we will explore the perspectives for open charm production in nucleus- 
collisions from AGS to RHIC energies employing the HSD transport approach |TO|JT^ for 
the overall reaction dynamics using parametrizations for the elementary production channels 
including the charmed hadrons D, D, D*, D*, D^, D„ D*, D*^, J/^, ^(25), X2c from A^A^ and 
7tN collisions. The latter parametrizations are fitted to PYTHIA calculations above y/s 
= 10 GeV and extrapolated to the individual thresholds, while the absolute strength of the 



cross sections is fixed by the experimental data p7|- ^] similar to Ref. (Section 2). The 
production of open charm in central collisions of Au + Au at 25, 160 A-GeV and 21.3 A-TeV 
is studied in Section 3 with respect to transverse mass (m^) and rapidity distributions. We 
will first switch off elastic collisions of kaons and 0-mesons with pions in nucleus-nucleus 
collisions to allow for a more transparent comparison to the spectra from pp collisions (Sect. 
3.1 - 3.3). The modifications of the m^-spectra due to the latter elastic interactions will be 
discussed in Section 3.4 while an excitation function for various mesons in central Au + Au 
collisions is presented in Section 3.5. Section 4 concludes this study with a summary and 
discussion of open problems. 



II. ELEMENTARY CROSS SECTIONS FROM PN AND vriV COLLISIONS 

Before examining nucleus-nucleus collisions we have to specify the differential open charm 
cross sections from pN and vrA^ reactions that will enter the HSD approach. Contrary to 
light meson production in hadronic reactions the creation of a cc pair is due to a hard process 
and dominated by gluon-gluon fusion at high ^/s. Using MRS G (next to leading order) 



structure functions from the PDFLIB package for the gluon distribution of the proton, 
a bare charm quark mass = 1.5 GeV and kx = 1 GeV we obtain the cross sections for 
D, D, D*,D*, Ds, Ds, D*, D*, as a function of > 10 GeV from PYTHIA ||] as displayed 



in Fig. 1 (upper part). Corresponding results for vrA^ reactions are shown in the lower part. 
Since the individual lines are hard to distinguish, some general trends are pointed out: All 
cross sections indicate a common (smooth) energy dependence. The D*-mesons are created 
more abundantly than the D-mesons roughly by a factor of 3 due to the three different spin 
polarizations; the small mass difference between D- and D*-mesons of ~ 140 MeV plays 
almost no role for y/s > 10 GeV. On the other hand, an exchange of a light {u, d) quark by 
a strange (s) quark costs a factor of 3-4. Consequently, the cross sections of D^, Z)°, D~ 
and D*,D* are roughly comparable at high ^/s. These simple considerations specify the 
relative abundance of the open charm mesons. However, the absolute magnitude of the 
cross sections is not expected to match experimental data due to the perturbative nature of 
these calculations and rescaling factors K have to be introduced . 

In this spirit we fit the individual results from PYTHIA (multiplied by factors of 12 and 
7 for pN and vrA^, respectively) by an expression of the form. 



axis = ax 



[1 - Zy Z"^, (1) 

with Z = \J^/ \/s where denotes the threshold for the channel X in pN or vrA^ 

reactions. Note that close to threshold the production for Z)(c)-mesons is enhanced relative 
to D(c)-mesons since the c-quark can end up in Ac, Sc, S* baryons with lower threshold, 
while D(c)-mesons require the associated production with a -D(c) meson. These threshold 
phenomena are in close analogy to the strangeness sector, where the mesons with a s-quark 
are produced close to threshold essentially together with hyperons (A, S), whereas antikaons 
require the associated production with a kaon. 

The formula (|l|) ensures the proper thresholds by construction while the exponents a 
and (3 describe the rise at threshold and the asymptotic behaviour, respectively. In order to 
properly 'normalize' the results from Fig. 1 we address to the experimental data from Refs. 



|37|-^ that have been extrapolated to full open charm cross sections by using the charge 



ratio's as given by PYTHIA. Furthermore, we have used a factor of 2 when extrapolating 
data for Xi? > to the full Feynman xp regime for pN collisions and a factor of 1.6 for vrA^ 
reactions ESI. The results of the fits are given in Tables 1 and 2 for the parameters ax,« 



and (3. We mention that the high value of the exponent a compared to related fits for p, uj 
or (f) production |]lOl indicates the different production mechanism for cc pairs compared to 



light quark pairs. 

The parametrized results from this extrapolation are displayed in Fig. 2 for pN (upper 
part) and vrA^ reactions (lower part) for the full charm cross section including all mesons 
as specified above with their individual thresholds. The solid lines in Fig. 2 represent the 
sum over all open charm mesons (within the parameters given in Tables 1 and 2) while the 
individual lines refer to the individual mesons that are somewhat hard to disentangle. As 
in Fig. 1 these cross sections group to 3 bunches at high ^/s where the upper bundle of 
lines corresponds to D*~ , D*^ and Z)*°, the middle bundle to D~ , _D° and the 
vector states with a strange quark D*, D*, while the lower bundle gives the cross section for 

It is interesting to compare these results with the cross sections for J/\E' (including Xc 
decay) and \E'' (^(25*)) which are displayed in Fig. 3 as a function of i/i together with 
the parametrizations (solid lines) for pN and vrA^ reactions (taken from Ref. [0,Q). Note 
that at ^/s = 20 GeV open charm is enhanced by about a factor of 50 relative to J/\I' and 
that this ratio increases with the available energy. Since the parametrization from Ref. ||50[] 
approaches some constant value at high ^/s contrary to the PYTHIA calculations (cf. Fig. 
1 of Ref. ||lT|) we have fitted the total cross section by the function 

axis) = bxil- YT 0(v^ - v/i^) (2) 



with Y = mxl\fs and a = 10, while ^/sq denotes the threshold in vacuum. Again the 
parameter (3 governs the high energy rise of the cross section which for /5 ~ 1 is now in line 
with the PYTHIA calculations specified above. Our fits give = 96 nb, b^^ = 64 nb, b^i 



= 20 nb; the results for J/'^ (including the Xc decay) are shown in the upper part of Fig. 3 
in terms of the dashed line. 

The cross sections (|T]),(0) will be used in the transport calculations to be discussed 
below which, apart from the total cross sections, also need the differential distribution of 
the produced mesons in the transverse momentum px and the rapidity y (or Feynman xp) 
from each individual collision. We recall that xp = Pz/vT"'^ ~ '^Vzj \fs with p^ denoting the 
longitudinal momentum. For the differential distribution in and pt we use the ansatz, 

^ , (1 - Vf\V exp{-bpT), (3) 

2pT dxpapT 

with 7 ^ 4.5 and b ^ 3.0 GeV~^. With these parameters the differential transverse momen- 



tum distributions oi D / D mesons in pp (and vrA^) reactions at 250 GeV [^] may reasonably 



be described as shown in Fig. 4. The xp and pt distribution for charmonium production, 



furthermore, is taken from Ref. pT| 



We have to point out that our parametrizations for the differential and total cross sections 
for open charm (as well as charmonia) become questionable at low energy, but also at high 
energy. It is thus mandatory that they have to be controlled by experimental data from pp, 
pA and vrA^ reactions before reliable conclusions on open charm dynamics in nucleus-nucleus 
reactions can be drawn. 

For the interpretation of the results from nucleus-nucleus collisions (cf. Section 3) it is 
worth to compare to pp collisions at different energies, respectively. To this aim we display 
in Figs. 5-7 the differential multiplicities {2mT)~^dNx / dniT in the transverse mass 



niT = ^Pt + m\ (4) 

for all final pions, kaons, 0-mesons, D + D mesons and charmonia from pp reactions at 
^/s = 7.1 GeV, 17.3 GeV and 200 GeV, respectively. The pion spectra describe the sum 
of vr"'", 7r°, 7r~, the kaon spectra the sum of , K^, , K~ , the D-meson spectra the sum 
of all D, D* , Ds, D* and their antiparticles while the spectrum denoted by cc includes the 
J/^P, the Xc as well as the ^P', where the latter contribution starts at rriT ~ 3.7 GeV and 



becomes visible as a tiny kink in the mT-spectra. Here the open charm and charmonia 
results stem from the parametrizations specified above (including the decay Xc J /'^ + 1) 
while the spectra for pions, kaons and 0-mesons are from the LUND string model (as 



implemented in the HSD transport approach). For orientation we also show exponential 
spectra with slope parameters of 143 MeV, 176 MeV and 225 MeV, respectively, which 
describe the m^-spectra of pions rather well. The kaon spectra at all energies are down by 
a factor of ~ 3, the spectra by a factor of 9-10 relative to this line due to strangeness 
suppression in pp collisions. However, it is quite remarkable that the charmonia spectra 
fit well to this approximate my-scaling (within a factor of 2-3) at = 7.1, 17.3 and 200 
GeV, respectively. Furthermore, the spectrum of open charm is roughly compatible with 
m-r-scaling at = 17.3 and 200 GeV, while the D mesons are suppressed relative to the 
scaling by a factor ~ 30 close to threshold {^/s = 7.1 GeV). Whereas these results basically 
stem from our parametrizations at ^/s = 7.1 and 200 GeV, the spectra at ^/s = 17.3 are 
controlled by experimental data. Such an 'apparent' statistical production of mesons in 
elementary reactions has been advocated before by Becattini [Q. 

We point out that also the approximate m-r-scaling from Figs. 5-7 has to be controlled 
by explicit experimental measurements. Data in a limited rapidity range might lead to 
somewhat different results since the rapidity distributions of pions, kaons, (/)'s, D's, D*'s 
and charmonia differ substantially due to kinematical reasons, i.e. the width of the rapidity 
distribution decreases with increasing meson mass. 

III. NUCLEUS-NUCLEUS COLLISIONS 

Inspite of the inherent uncertainties pointed out above it is worthwhile to explore the 
dynamics of open charm mesons in relativistic nucleus-nucleus collisions. Experiments are 
planned at the SPS as well as at RHIC and might be even performed in the 20-30 



A-GeV region Here we will employ the HSD transport approach for the nucleus- 

nucleus dynamics that has been tested in detail for pp, pA and AA reactions from SIS to 



SPS energies |T0|,|55|J56|] and been used for the description of charmonium production and 



propagation as well P^J57|J5^ . 

We recall that (as in Refs. [ p7| , pS[ ) the charm degrees of freedom are treated perturba- 
tively and that initial hard processes (such as cc or Drell-Yan production from NN collisions) 
are 'precalculated' to achieve a scaling of the inclusive cross section with the number of pro- 
jectile and target nucleons as Ap x At- To implement this scaling we separate the production 
of the hard and soft processes: The space-time production vertices of the cc pairs are cal- 
culated in each transport run by neglecting the soft processes, i.e. the production of light 
quarks and assosiated mesons. The resulting number Ncoii{b) of these 'hard' collisions is 
shown for Au + Au at 160 A-GeV in Fig. 8 (full squares) as a function of impact parameter. 
The inclusive number of inelastic A^A^ collisions is given by the integral of Ncou over impact 
parameter 

27iJbN,ou{b) db 2 

^ = — ~ ^ ' 

O'inel.W-'^) 

which gives approximately A^, i.e. the experimental scaling for 'hard' processes. In (^ 
the mass number A = 197 for Au, while crinei.{V^){^ 34 mb) denotes the inelastic nucleon- 
nucleon cross section. The calculated Ncoii{b) compares well with the result from Glauber 
theory (solid line in Fig. 8), where the number of inelastic interactions in nucleus- nucleus 
collision A+B at impact parameter b = {b, 0, 0) is given as 



NABib) = AB I a,neiTAis)TBis-h)d^s, (6) 

where s = {s^, Sy, 0) is orthogonal to the 2;- (beam-) direction. In the integral (^ 

+00 

TA{b) = J p{VWT^) dz (7) 

—00 

is the profile function normalized to unity, while p{r) is the nuclear density taken of Woods- 
Saxon shape. 

Thus the scaling for initial hard processes is adequately realized in the transport ap- 
proach. We mention that this scaling prescription might no longer be valid at low and high 



energy due to modifications of the gluon structure functions during the heavy-ion reaction 
or related shadowing phenomena pT|. For our initial study, however, we discard such effects. 



Apart from primary hard A^A^ collisions the open charm mesons or charmonia may also be 
generated by secondary 'meson'- 'baryon' reactions. Here we include all secondary collisions 
of mesons with 'baryons' by assuming that the open charm cross section (from Section 2) 
only depends on the invariant energy ^/s and not on the explicit meson or baryon state. 
Furthermore, we take into account all interactions of 'formed' mesons - after a formation 
time of Tp = 0.8 fm/c (in their rest frame) - with baryons or diquarks, respectively. 

In the transport calculation we follow the motion of the charmonium pairs or produced 
D, Z)-mesons within the full background of strings/hadrons by propagating them as free par- 
ticles, i.e. neglecting in-medium potential^, but follow their coUisional history with baryons 
and mesons or quarks and diquarks. For reactions with diquarks we use the corresponding 
reaction cross section with baryons multiplied by a factor of 2/3. For collisions with quarks 
(antiquarks) we adopt half of the cross section for collisions with mesons. Whereas the 
latter concept is oriented at the additive quark model, this assumption still does not solve 
the problem since the cross sections of D-mesons or charmonia with baryons and various 
mesons (essentially vr, p and uj mesons) are not well known. Thus we will provide results 
with and without rescattering of open charm mesons. 

In order to study the effect of rescattering we tentatively adopt the following dissociation 
cross sections of charmonia with baryons independent on the energy: 

accB = 6 mb; ctj/^b = 4 mb; a^^B = 5 mb; (t^>b = 10 mb, (8) 



while a lifetime (in it's rest frame) of 0.4 fm/c is assumed for the pre-resonance cc pair 
The energy-dependent J/^'-meson cross sections for dissociation to DD are taken from the 
calculations of Haglin which on average lead to a similar J/^P comover suppression than 
the overall cross section of 3 mb adopted in Ref. . 



^Except for the case of in-medium mass shifts in Section 3.2 



On the other hand, the D / ID mesons are expected to have large cross sections with mesons 
or baryons due to the hght flavor content such that light meson (vr, p,uj,ri) exchanges should 
describe the dominantly elastic cross sections at low relative momenta. We here adopt the 
calculations from Ref. |T^ which predict elastic cross sections in the range of 10-20 mb 
for D, D* scattering with mesons dependent on the size of the formfactor employed. As a 
guideline we use a constant cross section of 10 mb for elastic scattering with mesons and also 
baryons, although the latter might be even higher for very low relative momenta. We neglect 
charm exchange reactions such as D~^N ttAc, Sc or it Ac —>■ DN in the present study, which 
will essentially modify the charm quark content of mesons relative to baryons. Furthermore, 
we discard a recreation of charmonia by channels such as + Z) ^ J/^ + vr, since at AGS 



and SPS energies these reactions are negligible [30|; on the other hand, at RHIC energies this 



charmonium formation might become essential However, the formation cross sections 
are not well known and the significance of these channels is discussed controversely in the 
present literature PS|PT|,|B3I . 



In the transport calculations to be discussed below we will focus on the relative yield 
of pions, kaons, 0-, D + Z)-mesons and charmonia, that will be analyzed in terms of global 
m^-spectra which are integrated over the whole rapidity range. In order to allow for a more 
direct comparison with the m^-spectra from pN collisions in Figs. 5-7 we will first discuss a 
more transparent situation and 'switch off' a couple of reaction channels; the results from the 
'full' calculations will be presented in Section 3.4. To disentangle various dynamical effects 
we thus first suppress the decay (j) KK, irp to allow for a direct evaluation of the 0-meson 
my-spectra at the end of the calculation. Furthermore, it is well known experimentally 
that the apparent slope of my-spectra for different hadrons varies almost linearly with the 
hadron rest mass due to a common collective flow velocity (3. In the transport calculations 
this collective flow results from elastic collisions between the hadrons in the expansion phase 
of the reaction 1 34,351 ('pion wind'). For our exploratory study we switch off the elastic 
collisions of kaons and 0-mesons in the expansion phase with pions. This then leads to 
mr-spectra with roughly the same slope for all hadrons and their relative abundance can be 



extracted in a simple way (see below). As mentioned before, the related changes of the rriT 
spectra due to elastic collisions will be addressed in Section 3.4. 



A. SPS energies 

We now turn to the results of the HSD transport calculations. In Fig. 9 we show the 
time evolution of c, c production (solid line) for a central (b = 1 fm) collision of Au + Au at 
160 A-GeV in comparison to s, s production (dashed line), where the c, c number is scaled in 
height to the s, s line for the initial 'hard' production by a factor of 1.5 x 10^. Both functions 
rise steeply within about 1 fm/c; whereas the solid line (c, c) stays practically constant 
the dashed line (s, s) increases smoothly due to secondary and ternary ss production by 
meson-baryon or meson-meson collisions ||5^ . This 'cooking' of strangeness in the expanding 
'fireball' leads to a moderate (~ 46 %) enhancement of strangeness whereas the secondary 
production of cc pairs (by meson-baryon collisions) proceeds early and is only ~ 9% (for 
the cross sections specified in Section 2), which might be neglected at SPS energies. In this 
respect charm quark pairs dominantly are created in the initial high density phase of the 
collision with energy densities even above 3 GeV/fm^ [^. The multiplicity of open charm 
mesons here is about 0.2, whereas the multiplicity of J/\E''s (including the decay of Xc) is 
only about 10~^. The fraction of charmonia dissociated by baryons and mesons is 70% for 
J/\I', 80 % for Xc and 90 % for \E'', which is comparable to the suppression calculated earlier 



in Ref. 57 



The effect of rescattering of D-mesons on baryons and mesons as well as charmonium 
interactions with hadrons is shown in Fig. 10 for a central collision of Au + Au with respect 
to the transverse mass spectra. Here the D, D* spectrum is fiattened out in transverse 
mass while the charmonium spectrum is roughly reduced by a factor 3-4 due to dissociation 
reactions; the kink from \E'' at ~ 3.7 GeV disappears due to the large \E'' dissociation. We 
point out that a drastic enhancement of the slope of the D-meson m^-spectra as advocated 



in Refs. [32,33| is not seen from our dynamical calculations for the cross sections adopted. 



Quite remarkably, the open charm spectra and charmonium spectra appear to scale well in 
transverse mass after including the secondary interactions with hadrons. 

The effect of final state interactions on the rapidity distribution of D-mesons is displayed 
in Fig. 11 which shows a slight broadening of the distribution due to the elastic scattering 
processes with baryons and mesons, respectively. 

It is interesting to have a look at the my-spectra for all mesons in analogy to Fig. 6 
(for pp reactions) to explore the effects of open charm and charmonium rescatterings. The 
calculated mT-spectrum for pions, kaons, 0-mesons, all D + D mesons and charmonia is 
given in Fig. 12 which can be characterized well by an exponential slope parameter of 176 
MeV (dashed line) for all mesons. This result comes about as follows when compared to 
Fig. 6: The D-mesons are created more abundantly than pions (relative to pp) in central 
collisions of Au + Au because the D-meson (and charmonium) yield scales with the number 
of hard collisions (cf. Fig. 8) while the pions roughly scale with the number of participants 
and may be reabsorped to some extent. The kaon (and 0) yield increases due to rescattering 
(cf. Fig. 9), however, the (p {ss) mesons do not match the mT-scaling in the HSD transport 
approach and stay down by a factor of about 3-4. The charmonium spectrum (relative to pp) 
is decreased by a factor ^ 3-4 due to dissociation processes as noted before. All these effects 
lead to the approximate mT-scaling without employing the assumption of a Quark-Gluon 
Plasma (QGP) formation and a common hadronization at some temperature of 160 - 180 
MeV. 

We mention that a roughly constant vr to J/\E' ratio from pp to central Ph + Ph collisions 



has been also pointed out in Refs. [§^,0 proposing a statistical hadronization scheme in 
all reactions. Furthermore, Gallmeister et al. have suggested in Ref. |^ that the open 
charm degrees of freedom might be described in a simple thermodynamical model for central 
collisions of Ph + Ph at the SPS using the same temperature for all mesons. The findings of 
these authors are supported here by the nonequilibrium transport calculations, that provide 
rather simple arguments for the phenomena pointed out before. 



B. Au + Au reactions at 25 A GeV 



In this Section we explore the perspectives of open charm measurements in nucleus- 
nucleus collisions at 25 A-GeV, which might be accessable at a possible future GSI facility 
||5^ . In this initial study we restrict to central collisions of Au + Au at 25 A-GeV, which are 
expected to provide the optimal conditions for open charm experiments and studies on the 
in-medium properties of D-mesons in analogy to the K~^,K~ experiments at the SIS. We 
step ahead as in Section 3.1. 

Fig. 13 shows the time evolution of c, c production (solid line) for a central (b = 1 fm) 
collision in comparison to s, s production (dashed line) where the number of c, c is scaled 
again in height to the s, s line for the initial 'hard' production (by a factor 1.5x10^). Both 
functions rise within a few fm/c which corresponds to the passage time of the (Lorentz 
contracted) nuclei. As in Fig. 9 the solid line (c, c) stays constant for later times while 
the dashed line {s, s) increases again due to secondary and ternary ss production channels. 
The relative enhancement of ss 'cooking' here amounts to roughly 65% whereas the relative 
contribution of cc pairs from secondary channels is ~ 7% for the cross sections specified in 
Section 2. We note, however, that the 65% enhancement of strangeness is insufficient to 



explain the abundancies at the AGS |721 from 4-11 A-GeV or the K/tt ratio at 40 
A GeV (at the SPS) without assuming any in-medium modifications of the kaons. For a 
detailed discussion we refer the reader to Refs. [p^p6| . 

It is apparent from Fig. 13 that the charm quark pairs are created in the initial high 
density phase of the collision, here with energy densities up to 2 GeV/fm^, which is above the 
critical energy density from lattice calculations for the formation of a QGP [0. However, 
the energy densities from the transport calculation correspond to nonequilibrium phase- 
space configurations at high baryon density, that should not be identified with the energy 
density extracted from lattice calculations (in equilibrium and for quark chemical potential 
/i, = 0). 

For a quantitative orientation we display in Fig. 14 the volume (in the nucleus-nucleus 



center-of-mass) with an energy density above 1 GeV/fm^ and 2 GeV/fm^ as a function of 
time for a central Au + Au collision at 25 A-GeV, where only interacting and produced 
hadrons have been counted as in Ref. pB[. It is important to note that the high energy 
density is essentially build up from 'strings', i.e. 'unformed' hadrons. This phase may be 



addressed as string matter (cf. Refs. [^,^) and expresses the notion that most of the 
hadrons appear in some form of 'continuum excitation'. The energy density including only 
'formed' hadrons (during the expansion of the system) stays below 1 GeV/fm^, i.e. below 
the energy density expected for a transition to the QGP. The absolute numbers in Fig. 14 
have to be compared to the volume of a Aw-nucleus in the moving frame which, for a Lorentz 
7-factor of 3.78, gives ~ 330 fm^. Thus also at 25 A-GeV the phase boundary to a QGP 
might be probed in a sizeable volume for time scales of a few fm/c. Contrary to central 
collisions at the SPS these volumes are characterized by a high net quark density; for such 
configurations we presently have no reliable guide from QCD lattice calculations. 

The multiplicity of open charm mesons at 25 A-GeV is about 6 ■ 10~^, whereas the 
multiplicity of J/\E''s (including the decay of Xc) is about 1.5 ■ 10~^. We mention that the 
fraction of charmonia dissociated by baryons and mesons is ~ 60% for J/^ 

The effect of rescattering of D-mesons on baryons and mesons is displayed in Fig. 15 (for 
a central collision of Au + Au) for the D-meson rapidity distribution, which shows now a 
substantial broadening due to scattering processes with baryons and mesons. The decrease 
of the D-meson rapidity distribution at midrapidity is almost a factor of 2. 

The mT^-spectra for all mesons in analogy to Fig. 5 (for pp reactions) are presented in 
Fig. 16. The calculated mT-spectrum for pions, kaons and 0-mesons can be characterized 
by an exponential slope parameter of 143 MeV (dashed line). Again the kaon (and 0) yield 
is increased (relative to pp times the number of hard collisions Ncoii) due to rescattering (cf. 
Fig. 13), but the (ss) mesons stay down by a factor of 3-4. The charmonium spectrum 
(relative to pp) is decreased by a factor ^2.5 due to dissociation as noted before and 
approximately fulfills the global mr-scaling. The latter does not hold for D-mesons (open 
squares) which are suppressed dynamically in the threshold region by roughly one order of 



magnitude. 

We recall that a similar observation has been made for the m^-scaling of and K~ 
mesons close to threshold energies at the SIS |]73|, where the strange mesons have been 



suppressed relative to pions and ?7's. However, when adding to the K'^ mass the A — mass 
difference of 177 MeV (due to the associated production mechanism in pp and vrA^ collisions), 
a remarkable m-r-scaling could be recovered again |j7^. It should be noted that the latter 



scaling is not due to a grand-canonical (or canonical) chemical equilibration, but simply due 
to a shift of the spectra induced by the kaon production mechanism. We have to stress, 
however, that all these observations on the charm sector are based on our extrapolations 
(Section 2) and might not hold experimentally. 

We now address the question, to what extent in-medium modifications of the Z)-mesons 
might be seen in the m-r-spectra for central Au + Au collisions at 25 A GeV. Contrary 
to open charm production and propagation in antiproton induced reactions on nuclei ||26|| , 



where the D-mesons show up with momenta of a couple of GeV/c relative to the nuclear 
matter rest frame, the D-mesons produced in central nucleus-nucleus collisions have only 
small momenta in the rest frame of the hadronic fireball. This is of particular relevance 
for experimental studies of hadron self energies, since the latter are generally momentum 
dependent and most pronounced for low momenta. 

The modifications of the D-meson spectral functions in the medium at present cannot 
be reliably calculated nor extracted (in the low density limit) from experimental scattering 



data via a dispersion analysis (cf. Ref. ||75[ for the K^K problem). For our initial study 



we thus discard all momentum dependence of the D-meson self energies and also neglect a 



broadening of their spectral functions due to interactions in the medium |]7^. As a guide 
we employ the QCD sum rule calculations from Ref. and implement a mass shift of the 
form 

AniDip) = an— (9) 
Po 

with ~ —50 MeV, where po denotes the nuclear matter density and p the actual baryon 



density at the D-meson creation point. In principle one might expect different mass shifts 
of D and D mesons in the medium due to an opposite sign of the vector interaction [EB 



However, since the c, c pairs are created in the early high density phase of the collision 
(cf. Fig. 13) the vector interaction is expected to vanish here (cf. Refs. |]7^ , [77|] ) and only 
the scalar attraction to survive. This will lead to similar mass shifts for D and D mesons 
as anticipated in Eq. (H). Since densities up to 8po can be achieved in central Au + Au 
collisions at 25 A GeV, the D-meson mass shifts may reach up to -400 MeV. Such mass 
shifts have a dramatic effect on the production cross sections in pN collisions and secondary 
meson-baryon reactions when incorporating them in the production thresholds (cf. Tables 
1 and 2). 

Our calculations with the mass shift (crosses in Fig. 16) give an enhancement of the 
D-meson yield by about a factor of 7 relative to the bare-mass case (open squares). The 
slope of the spectra is modified only slightly relative to the bare mass case as can be seen 
from Fig. 16 for the resulting m-p spectrum. Somewhat surprisingly, an approximate itlt- 
scaling with all other mesons is regained in this case. We have to point out again that the 
results on open charm and charmonia in Fig. 16 essentially depend on our extrapolations in 
Section 2 and the assumed self energies (^), which are not controlled by data. On the other 
hand. Fig. 16 should be helpful in guiding the experimental analysis. 

C. Central collisions of Au + Au at ^/s = 200 GeV 

Apart from the low energy (threshold) regime, explored in Section 3.2, we also present 
predictions for central Au + Au collisions at RHIC energies of y/s = 200 GeV or 21.3 A-TeV, 
which will be investigated in the near future. For an overview of the predictions performed 
within the HSD approach we refer the reader to Ref. |T^. Here we extend our calculations 



to open charm mesons and provide transverse mass spectra for all mesons in analogy to 
Figs. 12 and 16. 

The multiplicity of open charm pairs for Au + Au at b=l fm and 21.3 A-TeV from the 



HSD approach is about 16, whereas the multiphcity of final J/\[^'s (including the decay of 
Xc) is only ~ 8 ■ 10~^ since in this case the dissociation by baryons and mesons is ~ 90% for 
J/\l' |T^. We mention that the production of charm pairs proceeds within less than 0.8 fm/c 
(in the nucleus- nucleus cms) and that the amount of cc-pairs from secondary meson-baryon 
reactions is approximately 11 % for the cross sections specified in Section 2. 

The rrtT-spectra for all pions, kaons and antikaons, 0-mesons, D and D mesons as well 
as charmonium states are presented in Fig. 17 in analogy to Fig. 7 (for pp reactions). The 
calculated rriT spectrum for pions, kaons, 0-mesons can be characterized by an exponential 
slope parameter of 225 MeV (dashed line). Similar to the lower bombarding energies in 
Sections 3.1 and 3.2 the kaon (and (p) yield is increased due to rescattering in the hot 
and initially high density mesonic fireball. The charmonium spectrum (relative to pp) is 
decreased by a factor ^10 due to dissociation as noted before and no longer fulfills the global 
mr-scaling. The D-mesons (open squares) are somewhat enhanced relative to the mT-scaling 
which, however, might be an artefact of the parametrizations in Section 2. On the other 
hand, the relative suppression of charmonia could be compensated by DD fiavor exchange 
reactions to J/\E' + vr etc., i.e. the inverse channels responsible for charmonium dissociation 



in interactions with mesons ||3T|. However, as mentioned before, these cross sections are 
presently not sufficiently known such that a final answer on the relative importance of these 
subsequent charmonium production channels has to wait for future. 



D. Collective acceleration of mesons in the expansion phase 

As mentioned above, the m^-spectra from Figs. 12, 16 and 17 provide a global view 
on the effect of chemical (inelastic) reactions, but should not be compared with experiment 
directly since the rescattering of kaons and mesons with pions has been switched off. 
In order to demonstrate the effect of the elastic scatterings, that have been discarded in 
Sections 3.1 - 3.3, we present in Fig. 18 the results for all mj- spectra from the transport 
calculations for Au + Au at 25 A-GeV, 160 A-GeV and 21.3 A-TeV now including the elastic 



rescatterings. Here the slopes of the pions are shghtly decreased whereas the slopes for kaons 
and 0-mesons increase at all bombarding energies. The spectra for the open charm mesons 
as well as charmonia do not change within the numerical accuracy since their rescattering 
with baryons and mesons had already been taken into account in the previous calculations 
in Sections 3.1 - 3.3. We note explicitly, that the high slope parameter for 0-mesons of ~ 
300 MeV seen experimentally at midrapidity in central collisions of Pb + Pb at the SPS 
by NA49 ||78| is not reproduced within the HSD calculations due to the weak coupling of 
0-mesons to non-strange hadrons. If this phenomenon is related to an early acceleration of 
strange quarks and antiquarks in a QGP phase or due to unexpected large rescattering cross 
sections is presently unclear. 

For comparison we discuss a scenario where the D, D mesons and charmonia emerge 



from a hadronizing QGP at rather low hadron density as advocated in In the latter 
scenario the collective expansion with velocity v± results from the pressure in the QGP phase 
and all hadrons freeze out at the same quark chemical potentials /iq,/is;/^c and temperature 
Ttherm- The slopcs of opeu charm mesons and charmonia then are expected to change 



as Tx ~ Ttherm + mx/2v'j_ for low momenta and as Tx ~ Ttherm y {'^ + v±)/i^ — v±) for 
momenta p± ^ rrix, where v± = (3 is a. collective velocity in the range 0.4 < f ^ < 0.6. Thus 
D, D mesons and charmonia {J/"^, '^') should show effective slopes larger than 0.35 GeV at 
SPS and RHIC energies. This conjecture might be tested soon experimentally and prove or 
disprove the moderate slopes of open charm mesons as predicted within the hadron-string- 
dynamics (HSD) approach. 



E. Excitation functions of mesons in central collisions 

In order to provide a more complete overview on meson production we show in Fig. 19 the 
excitation function of open charm mesons in central Au + Au collisions from AGS to RHIC 
energies without employing any self energies for these mesons. The i)-mesons with a c are 
produced more frequently at low energies due to the associated production with A^, Sc, 



similar to the kaon case (cf. lower part). At roughly 15 A-GeV the cross sections for open 
charm and charmonia are similar, while the ratio of open charm to charmonium bound states 
increases rapidly with energy. This behaviour is quite similar to the excitation functions in 
the strangeness sector when comparing , K~ and 0-mesons. Since the excitation function 
for open charm drops very fast with decreasing bombarding energy, experiments around 20 
A-GeV will be a challenging task since the multiplicity of the other mesons is higher by 
orders of magnitude. On the other hand, the perspectives for open charm measurements 
at RHIC appear promising since about 16 cc (or DD) pairs should be created in central 
Au + Au collisions according to our calculations. 

We mention that the excitation functions for the pions, kaons, eta's and (/)-mesons have 
been taken from Ref. |[T^, while the multiplicities for J/\l/ have been recalculated using the 



novel comover absorption cross sections from Section 2 and Ref. as well the parametriza- 
tion instead of the Schuler fit [^]. Since the numbers up to 500 A-GeV are compatible 
within 30% we do not discuss these differences in more detail. The higher J/\l/ multiplicity 
at RHIC energies is a direct consequence of the cross section (1^). 

IV. SUMMARY 

In this work we have calculated excitation functions for open charm mesons in central 
Au + Au reactions from AGS to RHIC energies within the HSD transport approach. The 
'input' open charm cross sections from pp and t^N reactions have been fitted to results from 
PYTHIA and scaled in magnitude to the available experimental data. In order to study 
the relative changes from central Au + Au to pp collisions, we have first switched off elastic 
final state interactions of kaons and 0-mesons with pions in order to suppress their common 
acceleration in the 'pion wind' during the expansion phase. Within the parametrizations 
and results from the LUND string model |B2[ - which is incorporated in the HSD approach 



- we find an mj-scaling for pions, kaons, D-mesons and J/\I^ in central collisions of Au + Au 
at the SPS (with an apparent slope of 176 MeV) without employing the assumption of a 



Quark-Gluon Plasma (QGP) and common freeze-out properties. We have shown that this 
result is essentially due to an approximate m^-scaling in pp collisions at -^i = 17.3 GeV and 
D,D and J/^ final state interactions. Furthermore, final state elastic scatterings change 
this conclusion to a moderate extent since the relative meson abundancies are not altered 
anymore and their spectra only get modified due to a common collective acceleration. 

At bombarding energies of 25 A- GeV a suppression of D-mesons by a factor of ~ 10 
relative to a global mT-scaling with a slope of 143 MeV is expected if no D-meson self 
energies are accounted for. On the other hand, attractive mass shifts of -50 MeV at po ^ 
when extrapolated linearly in the baryon density - lead to an enhancement of open charm 
mesons by about a factor of 7 such that an approximate mT-scaling for all mesons (cf. Fig. 
16) is regained. 

At RHIC energies of y/s = 200 GeV or 21.3 A-TeV the global mr-scaling is expected to 
hold also within a factor of 2-3 except for the charmonium states which - within the HSD 
transport approach - are dissociated by baryons and 'late comovers' to ~ 90%. On the 
other hand, the inverse reaction channels D + D ^ J/\E'+ meson etc. might lead again to 
charmonium enhancement as suggested in Ref. [0. Here we leave this question open for 
future analysis. 

However, as pointed out throughout this work, the elementary cross sections for open 
charm and charmonia in pp and vrA^ reactions have to be measured in the relevant kinematical 
regimes before reliable conclusions can be drawn in the nucleus-nucleus case. Experimental 
data in the 20 - 30 A- GeV with light and heavy systems will have to clarify, furthermore, if 
the quasi-particle picture of open charm mesons at high baryon density is applicable at all. 
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TABLES 



TABLE I. The parameters ax, a and (3 for pN reactions 



pN 


Meson 


[GeV] 


ax [mb] 


a 


f3 




5.605 


0.523 


4.92 


1.36 




5.069 


0.496 


4.96 


1.36 


D+ 


5.609 


0.469 


4.76 


1.40 


D- 


5.073 


0.363 


4.94 


1.44 




5.889 


1.775 


4.90 


1.34 




5.230 


1.275 


4.56 


1.42 


D+* 


5.896 


1.514 


4.64 


1.40 


D-* 


5.233 


1.384 


5.20 


1.36 


Dt 


5.813 


0.171 


5.12 


1.34 


D7 


5.373 


0.102 


5.58 


1.42 


Dt* 


6.101 


0.496 


4.88 


1.38 


D-* 


5.516 


0.283 


5.50 


1.46 



TABLE II. The parameters ax, a and /? for ttN reactions 



ttN 


Meson 


^ [GeV] 

V ^ L J 


ax [mbl 


a 






4.667 


0.273 


2.86 


1.28 




4.150 


0.247 


3.80 


1.26 


D+ 


4.671 


0.255 


3.22 


1.28 


D- 


4.154 


0.286 


3.50 


1.22 


D°* 


4.951 


1.076 


3.14 


1.22 


5°* 


4.292 


0.774 


3.80 


1.26 


D+* 


4.955 


0.719 


2.86 


1.32 


D-* 


4.296 


0.839 


3.40 


1.24 


Dt 


4.875 


0.0932 


3.62 


1.22 


DJ 


4.435 


0.0545 


3.70 


1.34 


D+* 


5.162 


0.284 


3.42 


1.24 


D-* 


4.578 


0.163 


3.64 


1.34 
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FIG. 1. The cross section for open charm mesons from PYTHIA [36| for pp (upper part) and 
vrA^ reactions (lower part) using MRS G structure functions, nic = 1.5 GeV and /ct = 1 GeV, 
respectively. The upper solid lines denote the sum over sdl D + D mesons. 




FIG. 2. The cross section for open charm mesons in the parametrization (1) using the parame- 
ters from Tables 1 and 2 in comparison to the experimental data from Refs. |^^- for pp (upper 
part) and ttN reactions (lower part). The upper solid lines denote the sum over all D + D mesons. 
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FIG. 3. The cross section for J/*!' and mesons in the parametrizations from Ref. |^] (sohd 
hnes) in comparison to the experimental data for pN (upper part) and vrA^ reactions (lower part). 
The J/^ cross sections include the decay from Xc mesons. The dashed line in the upper part shows 
the J/"^ cross section for the parametrization (|2|). 
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FIG. 4. The differential cross section for D/D mesons in transverse momentum (squared) for 
pp reactions at 250 GeV within the parametrisation (P) (sohd hne) in comparison to the data from 
Ref. [|5|. 
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FIG. 5. The transverse mass specta from pp collisions at Tiab = 25 GeV for pions (full squares), 
kaons (open triangles), and (/)-mesons (full rhombes) from the LUND string model [^] as imple- 
mented in HSD. The D + D meson (open squares) and charmonium (full dots) spectra - including 
the decay Xc ^ J /'^ + 1 ~ result from the parametrizations specified in Section 2. The dashed line 
shows an exponential with slope parameter Eq = 0.143 GeV. 
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FIG. 6. Same as Fig. 5 for pp reactions at Tiab = 160 GeV. The dashed hne shows an exponential 
with slope parameter £"0 = 0.176 GeV. 
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FIG. 7. Same as Fig. 5 for pp reactfons at ^/s = 200 GeV. The dashed hne shows an exponential 
with slope parameter Eq = 0.225 GeV. 
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FIG. 8. The number of hard cohisions Ncou as a function of impact parameter h in the HSD 
approach (fuh squares) for Au + Au at 160 A-GeV (see text) in comparison to the number of 
colhsions in the Glauber approach (sohd hne). 
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FIG. 9. Time evolution for the production of s, s (dashed line) and cc quarks (solid line, multi- 
plied by a factor of 1.5 x 10^) in the HSD approach for a central Au + Au reaction at 160 A-GeV. 
The c, c numbers have been scaled to the initial hard scattering processes. 
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FIG. 10. The transverse mass spectra of D + D mesons and J/^,^' mesons in the HSD 
approach for a central Au + Au cohision at 160 A- GeV. The open symbols denote the spectra 
without rescattering and reabsorption while the full symbols include the final state interactions. 
The thin dashed line shows an exponential with slope parameter Eq = 0.176 GeV. 
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FIG. 11. The rapidity distribution oi D + D mesons in the HSD approach for a central Au + Au 
coUision at 160 A-GeV. The dashed hne denotes the spectrum without rescattering and reabsorption 
while the solid line includes the final state interactions of D-mesons with hadrons. 
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FIG. 12. The transverse mass spectra of pions (full squares), kaons (open triangles), (/)-mesons 
(full rhombes), D + D mesons (open squares) and J/^', mesons (full dots) in the HSD approach 
for a central Au + Au collision at 160 A- GeV. The thin dashed line shows an exponential with 
slope parameter Eq = 0.176 GeV. Note that final state elastic scattering of kaons and </>-mesons 
with pions has been discarded in the calculations. 
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FIG. 13. Time evolution for the production of s,s (dashed hue) and cc quarks (sohd hue) in 
the HSD approach for a central Au + Au reaction at 25 A-GeV. The c, c numbers have been scaled 
to the initial hard scattering processes by a factor of 1.5 x 10^. 
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FIG. 14. Time evolution for the volume with energy density e > 1 GeV/fm^ (dashed line) and 
> 2 GeV/fm'^ (solid line) in the HSD approach for a central Au + Au reaction at 25 A- GeV. 
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FIG. 15. The rapidity distribution oi D + D mesons in the HSD approach for a central Au + Au 
colUsion at 25 A-GeV. The dashed hne denotes the spectrum without rescattering while the solid 
line includes the final state interactions of D-mesons with hadrons. 
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FIG. 16. The transverse mass spectra of pions (full squares), kaons (open triangles), (/)-mesons 
(full rliombes),D + D mesons (open squares) and J/^', mesons (full dots) in the HSD approach 
for a central Au + Au collision at 25 A-GeV without including self energies for the mesons. The 
crosses stand for the L)-meson spectra when including an attractive mass shift according to (|^. 
The thin dashed line shows an exponential with slope parameter Eq = 0.143 GeV. Note that final 
state elastic scattering of kaons and (/)-mesons with pions has been discarded in the calculations. 
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FIG. 17. The transverse mass spectra of pions (full squares), kaons (open triangles), (/>-mesons 
(full rhombes),D + D mesons (open squares) and J/^,^' mesons (full dots) in the HSD approach 
for a central Au + Au collision at 21.3 A-TeV. The thin dashed line shows an exponential with 
slope parameter Eq = 0.225 GeV. Note that final state elastic scattering of kaons and 0-mesons 
with pions has been discarded in the calculations. 
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FIG. 18. The transverse mass spectra of pions (full squares), kaons (open triangles), (/>-mesons 
(full rhombes),D + D mesons (open squares) and J/^', mesons (full dots) in the HSD approach 
for central Au + Au collisions at 25 A-GeV, 160 A-GeV and 21.3 A-TeV. The spectra at 160 A-GeV 
and 21.3 A-TeV have been multiplied by a factor of 10. The thin dashed lines show exponentials 
with slope parameters Eq = 0.143, 0.176 and 0.225 GeV, respectively. The deviations from the 
thin lines essentially are due to final state elastic scattering of kaons and (/)-mesons with pions, 
which accelerate the heavier mesons (see text). 
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FIG. 19. The multiplicity of D, D and J/^-mesons (upper part) for central collisions of Au+Au 
in the HSD approach including elastic and inelastic reactions, but no in-medium modifications of 
their spectral functions. The multiplicities for ,7], K'^ , and <j) (in the lower part) have been 



taken from Ref. |14] while the lines for D{c), D{c) and J/^ are the same as in the upper part. 



